Phase diagram as a function of temperature and magnetic field for magnetic 

semiconductors. 
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Using an extension of the Nagaev model of phase separation (E. L. Nagaev, and A. I. 
Podel'shchikov, Sov. Phys. JETP, 71 (1990) 1108), we calculate the phase diagram for degen- 
erate antiferromagnetic semiconductors in the T — H plane for different current carrier densities. 
Both, wide-band semiconductors and "double-exchange" materials, are investigated. 
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I. INTRODUCTION 



Degenerate antiferromagnetic semiconductors are ob- 
tained by strongly doping antiferromagnetic semiconduc- 
tors (e.g., europium chalcogenides or lanthanium man- 
ganites). Over a concentration range of doping impuri- 
ties, the ground state of these compounds will be a mix- 
ture of antiferromagnetic (AF) and ferromagnetic (FM) 
phases ^ |^. Starting from a pure compound and 
doping, we find that the ground state at T — H — 
is AF and non-conducting up to some conduction elec- 
tron concentration driven by impurities, Ud, from which 
the homogeneous state turns out to be unstable against 
phase separation. The ground state becomes then in- 
homogeneous with a simply connected AF phase and a 
multiply connected FM phase. The exact geometry of 
the multiply connected phase is of fractal nature. On 
increasing doping, a geometric transition takes place in 
which the topology of the sample changes due to perco- 
lation of the FM phase. Now, the ground state corre- 
sponds to a multiply connected AF phase and a simply 
connected FM one. The doping concentration at which 
this transition occurs is denoted by nx- At this point 
an important change in the conductivity is expected, the 
material becomes a conductor. On a further increase 
of doping, a certain concentration, n^, exits at which 
this phase-separated state starts to be unstable and the 
sample becomes again homogeneous but now FM. The 
physical reason for having a phase-separated state in an 
antiferromagnetic semiconductor at T = = is the 
dependence of the energy of the charge carriers that ap- 
pear on doping on the magnetic order of the lattice. The 
charge carrier energy is lower if they move in a FM re- 
gion. Therefore, by interaction with the spin system in 
the lattice, they are able to change its magnetic order, 
creating FM micro-regions and become trapped in them. 
In a degenerate magnetic semiconductor this is a cooper- 
ative phenomenon: a number of electrons are self-trapped 



in the same FM micro-region diminishing the energy per 
electron necessary to create it. When the carrier concen- 
tration is high enough, the phase-separated state turns 
out to be stable. For T, H ^ the free energy has differ- 
ent expressions for the AF and the FM parts. Because of 
this, a change in the temperature or the magnetic field 
produces a variation of the relative volume of the AF and 
FM parts. Therefore, a temperature- or magnetic field- 
induced percolation can occur, resulting in a complicated 
phase diagram. 

In this article, we calculate the phase diagram in the 
T — H plane for magnetic semiconductors within the dop- 
ing range < n < rtx- The phase diagram is expected 
to be divided into three different stability regions, those 
corresponding to: insulating phase-separated state, con- 
ductive phase-separated state and homogeneous state. 
Our calculations are based in a variational principle for 
the free energy of the system developed by Nagaev and 



Podel'shchikov in the reference | 
approach of the references ||l|, ||, 
temperature. 



that generalizes the 
to the case of finite 
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II. CALCULATION OF THE FREE ENERGY 

In this section, we briefly review the variational 
method used. More detailed calculations can be found in 
references |l|, |, |, |. The microscopic description of the 
sample is provided by the Hamiltonian of the generalized 
Vonsovsky s-d model. The main parameters in this model 
are W = 2zt, the carriers band-width, AS the exchange 
energy between conduction electrons and magnetic ions 
(s-d exchange energy) , and zIS'^ the exchange energy be- 
tween magnetic ions (d-d exchange energy). / < is the 
exchange integral between first-nearest neighbors mag- 
netic ions. The smallest parameter is the d-d exchange 
energy. We differentiate two possibilities depending on 
the relative value of W and AS. In the case W >> AS, 
we have a wide-band semiconductor (e.g. EuTe). In 
the opposite case W « AS, we talk about a "double- 
exchange" material (e.g. lanthanium manganites). For 
the sake of definiteness the sign of A is assumed to be 
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positive. 

The free energy is expresed as a function of two vari- 
ational paramaters: the ratio x of the volumes of the 
AF and FM phases, x — y^, and R, the radius of the 
spheres which form the multiply-connected phase. All of 
the electrons are assumed to be found in the FM part 
of the phase-separated state Although this is not 
strictly true at T 0, it is a good approximation be- 
cause the number of electrons in the AF part is propor- 
tional to eT^, where U W ii W » AS or U AS 
if << AS, and in both cases we have U » T. The 
trial free energy is given by the expression: 

F = El^in + Esur +Ec+ F„,ag (1) 

where -Ekin is the standard bulk kinetic energy of the con- 
duction electrons, -Esur is the surface electron energy, i.e. 
the correction to -Ekin that appears due to the quantiza- 
tion of the electron motion in regions of finite dimensions 



. In view of the degeneracy of the electron gas, the con- 
tribution of the thermal excitations to its free energy can 
be neglected. Ec is the electrostatic energy due to the in- 
homogenous electronic density and Emag is the magnetic 
free energy calculated in the mean field approximation. 
Due to the absence of conduction electrons on the AF 
part of the phase-separated state, this trial free energy is 
vahd for both W » AS and W « AS. 

The expressions for these quantities are: 

3 2 
Evin = -^i{n)n{l + x)^ 
5 
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(2a; -I- 3 — 3 (1 -f x) ^ if FM spheres in an AF matrix, 

a; ^2 + 3a; — 3x^ (1 + x) x^^ if AF spheres in a FM matrix. 



where fi (n) = t (6n^nj ^ and /3 = 3 if the FM part is 
multiply connected, ot (3 = 3x if it is simply connected. 
The Coulomb energy per unit cell is found by separation 
of the crystal into Wigner cells. This method provides a 
good approximation to the electrostatic energy for small 
volumes of the alien phase, but it supposes to admit that 
the crystal is isotropic and homogeneous regarding the 
spatial distribution of these spheres. So we cannot de- 
scribe the fractal nature of the boundary that separes 
both phases. 



The free energy of the magnetic system is calculated 
using the mean field approximation. It is assumed that 
in the FM part conduction electrons are completely po- 
larized. Then: 



Ernag = — — E^ag(Er, T) + -_F„,ag(i^ + H,, T) (4) 

i I X J- I X 

where: 



Emag {H, T) 



^Sl + ljSf-Tln 
-\JSl-T\n 



if iJ < 2 J Si, 
if iJ > 2 J Si. 



(5) 



and the values Si, S2 verify the following equations: 
Si = SBs 

S. = SBs (^^1^ 



(6) 



where Bs is the Brillouin function. Moreover, we have 



defined He = "^■^"^^+^) as the mean field created by the 
conduction electrons on the localized d-spin system. The 
preceding energies where obtained taking into account 
that in the case of small fields, H < 2 J Si, the localized 
d-spins system is in a non-collinear AF (NCAF) state, 
being 5*1 the mean d-spin value in each of the sub-lattices. 
In the case of high fields, H > 2 J Si, the sample goes to 
a non-saturated FM (NSFM) state, being S2 the mean 
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value of the localized d-spin per site . 

III. RESULTS. 

The stationary state of the system is determined from 
the condition that the total free energy of the system, 
equation be a minimum with respect to the varia- 
tional parameters x and R. The parameter R only ap- 
pears in the surface and Coulomb energies. Minimizing 
the sum Q = E^ur + Ec with respect to R for fixed x 
leads to the following expression: 

Qmi„ = a7/?(|)^ (7) 
where the following definitions are used: 




After substitution of equation (|7|) in equation min- 
imization with respect to x is carried out numerically. 

For n < TiT three different phases are expected: insu- 
lating phase-separated, conducting phase-separated and 
homogeneous states. We calculate the regions of abso- 
lute stability of each phase in the plane T — H for dif- 
ferent conduction electron concentrations in the sample. 
The boundary lines between phases are calculated nu- 
merically. Below we present these diagrams for different 
concentrations, all below nx, for both wide-band Q and 
"double-exchange" materials. 

Our results indicate that for both kinds of materials, 
all of the phase transitions between homogeneous and 
phase-separated states are first-order. Preliminary cal- 
culations including the possibility of spin waves confirm 
this result . The transition between insulator and con- 
ductive phase-separated states is found close to second- 
order. The jump in the parameter x at the transition is 
about 20% at T = 0, growing with the transition temper- 
ature. This small jump in the parameter does not permit 
elucidating whether in a real sample the transition is of 
first- or second-order. We have to stress that the bound- 
ary between the separate states is of a fractal nature, and 
our approximation, that considers it as homogeneous and 
isotropically distributed spheres, does not work near the 
transition point. In any case, due to the smallness of the 
jump in x our approximation provides accurate values of 
T and H at the transition point. 

A. W » AS 

Phase diagrams for wide-band semiconductors are 
shown in figures (|^), (|^) and (^) for conduction elec- 
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FIG. 1: Phase diagram for a wide-band AFS with carrier 
concentration n = 0.10 riT 
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FIG. 2: Phase diagram for a wide-band AFS with carrier 
concentration n = 0.50 nx 
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FIG. 3: Phase diagram for a wide-band AFS with carrier 
concentration n = 0.95 nx 



tron concentrations equal to n = 0.10, 0.50, 0.95 nx, re- 
spectively. The areas 1, 2, 3 denote regions of abso- 
lute stability for the homogeneous state, the conductive 
phase-separated state and the insulating phase-separated 
state, respectively. The parameters that characterize 
these materials are chosen to be the same as in refer- 
ence [Q (they correspond to rare-earth compounds like 
europium chalcogenides): S = |, \J\S — 10^^ eV (im- 
plying Tn = 5 K, and the field at which both sublattices 
collapse is 98.7 kOe), AS = 1 eV, = 20, a'^ = 4 • lO^^ 
cm"'^, and the electron effective mass is equal to free elec- 
tron mass, this implies W = A eV. The value obtained 
for nx = 1.10 • 10^" cm"'^. The value for nu, at which 
the conducting phase-separated state becomes unstable 
at T = ff = 0, is 1.79 • 10^° cm'^. 

We see that at low carrier concentrations only a tran- 
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sition from the insulating phase-separated state to the 
homogeneous state is possible (see figure (|l|)). The phys- 
ical explanation of the diagram (|l|) is the following. At 
fixed n, the effect of the increase of magnetic field and 
temperature is to decrease the value of x with respect to 
its value at T = H = 0. In the range of concentrations 
we are working, rid < n < ^t, the state of the sample 
at r = = is an insulating phase-separated state 
{x > 1). Regarding to the behavior of a; at T = iJ = 
with n, X — > 1 as n — > nx- Following equations (|^), 
the d-spins of the magnetic ions in the FM part of the 
phase-separated state experience a magnetic field which 
is the sum of the external field, H, and an additional 
field created by conduction electrons. He. If we set the 
external field equal to zero, He is large enough to es- 
tablish the FM order in the multiply-connected part, i.e. 
He > 2 J Si. But as we increase external field, the value of 
X decreases very quickly and it is reached a point in which 
H + ^(1 + x) = 2JSi. On further increase of H , the 
magnetic field cannot maintain the multiply-connected 
part in a FM (NSFM) state and a transition to an ho- 
mogeneous AF (NCAF) state occurs. The transition line 
can be drawn in a simple approximate way. We start by 
looking for the value of x, xq, which makes the energy 
minimal at T = 0. For n = O.lOnx we found xq — 2.525. 
Because this value increases only slightly as we move 
along the transition line, this can be approximated by 
the expression H + ^{1 + xq) ^ 2JSi. 

At T = for n > 0.18 nx, the carrier concentration is 
high enough to permit that increasing the external field 
we remain in the insulating phase-separated state until 
X = 1. Therefore for concentrations n > 0.18 nx, three 
different phases appear. At low temperatures and low 
magnetic fields, we have a transition between both of 
the phase-separated states. This line correspond to the 
simultaneous percolation of the magnetic order and the 
conduction electron system. It is almost a second order 
transition. The values of x at both sides of the transition 
line are x — 1.09 and x = 0.92 at T = not varying 
with n, increasing as we move along the line towards 
higher temperatures, and being the growth smaller, the 
closer n to nx. To calculate the values of T,H,x^,X2 
(being 2:3, X2 the values of x at both sides of the transition 
line), we impose three conditions: the free energy of both 
phase-separated states must have a minimum and the 
free energy of both states must be equal. This allows 
us to calculate H,X3,X2 varying T along the line. Over 
the same range of temperatures but at higher fields, we 
find a transition between conductive phase-separated and 
an homogeneous state. This is a first order transition, 
similar to the melting of a solid. Again the same effect 
than in diagram (|^) appears once we are in the conductive 
phase-separated state, the value of n must be high enough 
to maintain the simply-connected in a FM state. If not, 
a first-order transition to an homogeneous state with a 
value of X approximately constant along the line occurs. 

In figure (^) , n is high enough to permit that increasing 
the magnetic field we remain in the conducting phase- 



separated state until a; = 0, so this effect is absent. In 
this case the transition line is determined in the same 
way as in the usual melting transition ||^. To calculate 
the value of r, H, X2 , we impose the conditions that the 
first two derivatives of the free energy of the conductive 
phase-separated state be equal to zero. At higher tem- 
peratures we find the transition line between the insulat- 
ing phase-separated and homogeneous states, that now 
corresponds to a NSFM order. It is important to notice 
that what we denote as phase 1 in the diagrams, cor- 
responds only to homogeneous magnetic states, without 
distinguishing whether they are of the NCAF or NSFM 
type. This is the reason why a kink in the transition line 
appear at T ~ 5 K. Really, it corresponds to a tri-critical 
point between insulating phase-separated, homogeneous 
NSFM states (at lower temperatures) and homogeneous 
NCAF (at higher temperatures). The transition line be- 
tween these two is not drawn. The transition line from 
the kink to the T-axis correspond to the magnetic tran- 
sition of the AF part of the insulating phase-separated 
state to a PM state. With n — > nx, this tri-critical 
point tends to the point T = T^, H — 0. Due to the ab- 
sence of electrons in the AF part of the phase-separated 
states the whole the sample has the same Neel temper- 
ature, as is confirmed by experiment. This part of the 
line is calculated in the same way that the transition line 
between conductive phase-separated state and homoge- 
neous states. 

The third diagram corresponds to n = 0.95 nx. We 
find two transition lines. One between both of the phase- 
separated states, close to second order, as in the pre- 
ceding discussion. We see that now the tricritical point 
between conductive phase-separated, insulating phase- 
separated and homogeneous states is absent. Therefore 
the percolation with temperature (this is, at H = 0) is 
also possible. The other line corresponds to the transition 
between conductive phase-separated and homogeneous 
states (now, NSFM all along the line). The kink was 
moved to higher temperatures. We see that the region of 
absolute stability of the insulating phase-separated state 
has been reduced, tending to dissappear as n goes to nx. 



B. W«AS 

Phase diagrams for "double-exchange" semiconductors 
are essentially the same. We also present three of them, 
in figures (||), (||), and (||), corresponding to conduction 
electron concentrations equal to n = 0.10, 0.50, 0.95 nx, 
respectively. The parameters that characterize "double- 
exchange" semiconductors are the following (typical for 
lanthanium manganitcs, for example) j9|: 5 = 2, \ J\S = 
2 ■ 10^2 cV (that implies Tn = 116 K, and the field 
at which the two sublattices collapse at T = is 3.46 
MOe), er = 5, a^^ = 4 • 10^^ cm~^, the effective electron 
mass equal to the free electron mass, that implies W — 4 
eV and A = A eV. With these parameters we obtain 
nx = 9.16 • 10^0 cm-3 and = 1.47 • lO^^ cm'^. 
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FIG. 4: Phase diagram for a "double exchange" AFS with 
carrier concentration n = 0.10 riT 
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FIG. 5: Phase diagram for a "double exchange" AFS with 
carrier concentration n = 0.50 jit 
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FIG. 6: Phase diagram for a "double exchange" AFS with 
carrier concentration n = 0.95 tit 



The behavior of the "double exchange" compounds is 
quantitatively different. The scenario discussed above 
for wide-band case remains valid, but a higher value of 
A implies a very strong magnetic field acting on the 
FM part of the separate states and also a higher value of 
m. It must be stressed that the values that we obtain for 
riT) are far from the true theoretical values. In particular, 
the values obtained depend on the physical parameters 
such as while we are dealing with a geometric tran- 
sition that cannot depend on these details. This is be- 
cause we lost the fractal nature of the surface that separes 
phases to simplify calculations. The quantum nature of 



the d-spin system forces the true percolation concentra- 
tion to be ripor = 0.16a~'^, the percolation concentration 
for the Bethe lattice (see, for example, indepen- 
dently of the physical processes involved. So, in both 
cases treated here, it is necessary to "map" rtx — > ?^pcr 
in order to compare with the experimental values. Be- 
cause HcS >> Tn we can neglect the terms with m ^ S 
in the partition function of the FM part. This makes 
the FM part temperature-independent, remaining in the 
same state as at T = 0. We see in figures (||) and (||), 
that the tri-critical point corresponding to the kink at 
T ~ 100 K, is now between conductive phase-separated, 
homogeneous NSFM (at lower temperatures) and homo- 
geneous NCAF (at higher temperatures), and appears 
before in temperature than the tri-critical point between 
conductive phase-separated, insulating phase-separated 
and homogeneous state (now, NCAF). As we see in fig- 
ure (^, percolation is possible varying the magnetic field 
in the same way as above, but percolation in temperature 
(i.e. at H ~ 0) is not, contrary to what was seen above. 
However, this result can be modified beyond a mean field 
approximation. It must also be stressed that the value 
of AS sets the scale of the H-axis and its high value 
makes the transition to the homogeneous states occur to 
unphysical values of the external field. 

A comparison with the experimental results is spe- 
cially relevant for the case of "double-exchange" com- 
pounds. It must be stressed that although all along this 
article we talk about the charge carriers as being con- 
duction electrons, the model is general enough to deal 
also with hole-doped compounds. An example of the lat- 
ter could be the canonical CMR lanthanium manganite 
Lai-xCaajMnOs with x < 0.16, which is extensively re- 
viewed in the reference However the present calcu- 
lation is more interesting for electron-doped compounds, 
such as Smi_a;Caa;Mn03 with x > 0.84, because the 
region of record negative magnetoresistance fall within 
the range of doping studied in this work. This com- 
pound is studied in the reference , where the observed 
CMR effects are explained as consequence of the metal- 
insulator transition described here. They also present 
the phase diagram in the T — H plane for th compound 
Smo.isCao.ssMnOs, which agrees qualitatively and quan- 
titatively with our figure (^. 
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